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rn . ABSTRACT 
O 

o We explore acceleration of ions in the Quark Nova (QN) scenario, where a 

neutron star experiences an explosive phase transition into a quark star (born in 
the propeller regime). In this picture, two cosmic ray components are isolated: 
one related to the randomized pulsar wind and the other to the propelled wind, 
both boosted by the ultra-relativistic Quark Nova shock. The latter component 
^ ■ acquires energies 10 15 eV < E < 10 18 eV while the former, boosted pulsar wind, 

achieves ultra-high energies E > 10 18 ' 6 eV. The composition is dominated by 
ions present in the pulsar wind in the energy range above 10 18 ' 6 eV, while at 
energies below 10 18 eV the propelled ejecta, consisting of the fall-back neutron 
star crust material from the explosion, is the dominant one. Added to these 
two components, the propeller injects relativistic particles with Lorentz factors 
Tprop. ~ 1 — 1000, later to be accelerated by galactic supernova shocks. The QN 
model appears to be able to account for the extragalactic cosmic rays above the 
ankle and to contribute a few percent of the galactic cosmic rays below the ankle. 
We predict few hundred ultra-high energy cosmic ray events above 10 19 eV for 
the Pierre Auger detector per distant QN, while some thousands are predicted 
for the proposed EUSO and OWL detectors. 



Subject headings: acceleration of particles — cosmic rays — elementary particles 
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1. Introduction 

Large efforts have been devoted to explore the origin of cosmic rays. Most puzzling are 
the observed ultra-high energy cosmic ray events (UHECRs) above the GZK-cutoff (Greisen 
1966; Zatsepin & Kuzmin 1966): protons lose energy drastically due to pion photo-production 
processes on the cosmic microwave background (CMB) at energies higher than 7 x 10 19 eV, 
and this limits the proton mean free path to some tens of megaparsecs (see Blasi (2003) and 
Bahcall & Waxman (2003) for more discussion on the GZK cut-off). Two distinctly different 
classes of models are commonly considered, the top-down and bottom-up scenarios. In the 
former ones, UHECRs are associated e.g. with the decay of some supermassive particles, 
whereas in the latter scenarios UHECRs are assumed to be accelerated by astrophysical 
objects (e.g. Sigl et al. (1994); Venkatesan et al. (1997); Biermann (1998), and for a latest 
review we refer the interested reader to Ostrowski (2002); Arons (2003); Sigl (2003); Olinto 
(2004)). A growing number of bottom-up explanations have been proposed, including active 
galactic nuclei, gamma ray bursts and neutron stars. However, at present there seems to 
be no clear association of UHECR events with any of these objects although by tuning the 
model parameters one enables the highest energy events to be accounted for. 

Identifying sources of the observed UHECR events remains uncertain and debated leav- 
ing room for speculation (Nagano & Watson (2000); Torres & Anchordoqui (2004); Stanev 
(2004)). As such we wish to consider in this paper the possibility that hypothetical quark 
novae (hereafter, QNe; Ouyed, Dey & Dey (2002)) contribute to the cosmic ray flux, es- 
pecially above 10 18 eV. In the QN explosion the core of a neutron star (NS) shrinks into 
the equilibrated quark object/star (QS). The overlaying crust material free-falls following 
the core contraction releasing enough energy to form an ultra-relativistic ejecta. The ultra- 
relativistic shock interaction with its surroundings environment (namely the randomised 
relativistic wind of the progenitor) leads to the first cosmic ray component in our model. 
The compact remnant acting as a magnetohydrodynamic propeller provides the second com- 
ponent, injecting ions to be boosted by the ultra-relativistic quark nova shock. Most of the 
propelled wind of relativistic particles is injected into interstellar space with low energies 
(from GeV to TeV), and these particles can be further on accelerated by galactic supernova 
shocks. Thus the propeller can contribute to cosmic ray flux also below the knee as an 
injector. It appears that the QN model can largely contribute for the extragalactic cosmic 
rays above the ankle and contribute a few percent of the galactic cosmic rays around and 
below the knee. 

The main assumption in this paper is the conversion of a highly magnetised, rapidly 
spinning NS to a QS. This phase conversion (or second explosion) - which remains to be 
confirmed - leads to unique conditions where two UHECR components are feasible. It is this 
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unique aspect of the model that makes this investigation worth pursuing. We note however 
that the approach presented here (including the acceleration mechanisms as related to the 
dynamics of the explosion) borrows heavily from what has already been presented in the 
literature for models of UHECRs in the context of isolated pulsars and binary coalescence 
(e.g., Venkatesan et al. (1997); Gallant & Achterberg (1999); Blasi (2000); Arons (2003) to 
cite only few). 

This paper is presented as follows: We start with a brief review of the concept of quark 
nova and its features (§2). In §3 we describe conditions under which the compact remnant 
is born in a propeller regime. In §4 we explain how cosmic rays are produced and isolate 
the two components (with energies 10 15 eV < E < 10 18 eV and E > 10 18 ' 6 eV). Here we 
describe the QN compact remnant as an injector of relativistic particles to the galaxy. In §5 
we predict UHECR events in future detectors and discuss multiple events and the isotropy 
of sources in our model. A discussion and a conclusion follow in §6. 

2. Quark Nova 

It has been suggested that the core of neutron stars may deconfine to a composition of up 
(u) and down (d) quarks during or shortly after some supernova explosions when the central 
density of the proto-neutron star is high enough to induce phase conversion (see e.g. Dai & 
Peng (1995); Xu, Zhang, & Qiao (2001)). It has also been speculated that when the density 
in the core increases further, a phase with strange quarks (s) becomes energetically favoured 
over the pure (u,d) phase and soon the entire star is contaminated and converted into this 
(u,d,s) phase. This is one of the scenarios introduced to convert an entire NS to a QS (e.g., 
Cheng & Dai (1996); Bombaci & Data (2000)). In the QN picture the (u,d,s) core is assumed 
to shrink to a corresponding stable quark object before the contamination is spread over the 
entire star. By physically separating from the overlaying material (hadronic envelope) the 
core drives the collapse (free-fall) of the left-out matter leading to both gravitational energy 
and phase transition energy release as high as Eq^ ~ 10 53 ergs (Ouyed, Dey & Dey (2002); 
Keranen, Ouyed & Jaikumar (2005)). 

2.1. Quark Nova ejecta and compact remnant 

The QN-ejecta consist mainly of heavy nuclei (the NS crust), protons, electrons and 
neutrons. The total amount of the ejecta has been estimated to be of the order of (0.001 — 
0.01)M Q with corresponding Lorentz factor Tq N = e Eq N / M C j CC c 2 ~ 10 — 100, where e is 
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the percentage of the QN energy transferee! to the QN shock front which we assume to 10%. 

The radius of the newly formed QS is given by R ~ -RnsIpns/p) 1 ^ 3 where the NS 
core density, p NS , and the QS density, p, scale as Pns/p — 0.1 — 0.2 (Ouyed, Dey & Dey 
2002; Keranen, Ouyed & Jaikumar 2005) - hereafter parameters with no subscripts refer 
to the quark star. The QS spins up during the phase transition due to contraction, with a 
rotational period P = P ns (-R/-Rns) 2 -Pns(Pns/p) 2 ^ 3 — -Pns/4- In the process, the magnetic 
field is amplified, B = B NS (R NS /R) 2 = 5 N s(p/Pns) 2/3 ~ 45ns. 



2.2. Fall-back material 

Fall-back of material onto a newly formed quark star may take place (similar to what 
has been suggested in the supernova case; Chevalier (1989)). In the early stages the accretion 
rate is given by 

m ~ 10 28 e/s ( — )( - ) 3/2 ( M ) 1/2 (1) 

where pg is the average density of fall-back matter (10 6 g/cm 3 , representing the crust mate- 
rial, and the matter below the neutron drip line densities). The hyper-Eddington accretion 
rate given above is understood by noting that (i) the fall-back matter in the initial phase of 
the explosion is not accreted onto the surface (but is propelled away before radiating), and 
(ii) the crust would not form in the early stages leaving the quark star bare (not subject to 
the Eddington limit since the bulk of the star is bound via strong interaction rather than 
gravity; see e.g. Alcock, Farhi, & Olinto (1986); Zhang, Xu k Qiao (2000), and Xu (2003) 
for a recent discussion on the matter). 



3. Propeller regime 

For the remainder of this paper, we consider a QS (NS) mass of M — 1.5M© (M NS ~ 
1.5M Q ), a radius R — 10 km (-Rns ~ 12.5 km), a surface magnetic field of 4 x 10 14 G 
(_B NS ~ 10 14 G). Since the fastest pulsar has a period of 1.56 ms (Backer et al. 1982), we 
use 2 ms as a representative value for the period P of the new born quark star (Pns ~ 8 
ms). Clearly these are unique constraints (millisecond, highly magnetised NS) favouring a 
scenario where the QS forms immediately following a SN explosion, or that where the NS 
has been spun up by accretion from a companion. 

The newly born quark star is defined by its three critical radii: the Keplerian "co- 
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rotation radius" 

R c = 27\nn(-^— ) 1 /3 ( JL ) 2/3 (2) 

the magnetospheric radius at which the ram pressure of the in-falling matter balances the 
magnetic pressure 

_ B 2 R 6 2/7 _ B 4/7 

^ m " ( 2mv / 2GM ) " ( 4 x 10" G } X 

I0 28 g/S 2/7 i? .12/7, 1-5 M 1/7 

(see, e.g. Frank, King, & Raine (1992)), and the light cylinder radius 

fl '<=S = 96km( 2^- (4) 

Given our fiducial values, the QS is born in the propeller regime (Shvartsman 1970; Illarionov 
& Sunyaev 1975), i.e. R c < R m < R\ c , where the infalling material may be accelerated in 
a wind that carries away angular momentum from the magnetosphere and hence from the 
QS itself (§4.2). Note that the propeller regime is only possible if the quark star is born 
with a magnetic field exceeding the critical value B c > 10 14 G, which corresponds to a 
parent neutron star initial magnetic field B^s, c — B c /A > 2.5 x 10 13 G. Here "c" stands for 
"critical" while the factor 4 (derived from our fiducial values) is the result of the amplification 
of the magnetic field following the quark-nova explosion. These critical values explain our 
assumption of a highly magnetised, rapidly spinning NS. 



3.1. Propeller lifetime 

The star loses rotational energy at a rate defined by the gravitational radiation losses 
(-Eg rav .), electromagnetic radiation losses (E em ) and the propeller's torque (E pTop .). That is, 

c\Q ■ E -+- E + E 

dt 2 m (5) 

with (Shapiro & Teukolsky 1983) 

-£ g rav. = - 5 = 6.8 x 10 47 erg/s x 
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where I is the moment of inertia and e the equatorial eccentricity; and with (Manchester & 
Taylor 1977) 

4P 2 P 6 ft 4 4fi , 

-E em = — — 2.6 x 10 46 erg/s x 

X (_ P~ )4( 4x 10 14 G )2( 10kn^ )6 - (7) 
The spin-down due to the propeller is (Menou et al. 1999) 

-E piop . = 2mc 2 = 1.8 x 10 49 erg/s (^^) ■ (8) 

The expression above assumes that the material flung away by the propeller effect has been 
accelerated to an angular speed corresponding to that of the star. We note that gravitational 
losses would be important in the spin-down for very short periods P < 2 ms (not our case) 
and propeller losses dominate over the electromagnetic dipole radiation losses given our 
accretion rates. 

The propeller regime lifetime can be obtained using Eq. (5) and Eq. (8). We find, 
assuming a constant accretion rate, 

t ~1Q3 S ( 1 )( 1 ° 28gS ' 1 )(^-^) (9) 

where P ; and Pf, initial and final periods, are in milliseconds. For a constant accretion rate, 
and taken the total amount of the QN ejecta, the lifetime of the propeller phase would not 
exceed a hundred seconds during which time the QS would have spun-down by no more 
than 30% thus remaining within milliseconds period. Following this phase, the spin-down is 
governed by the magnetic dipole radiation losses of the QS. 



4. Production of the two cosmic ray components 

4.1. Acceleration of the pulsar wind 

We assume that the pulsar wind is composed of electron-positron pairs and ions, as 
would be expected if pulsars operate as open-circuited systems with ions carrying the return 
current (Hoshino 1992). Some fraction of ions can be heavy like iron, originating from the 
surface of the pulsar. We also assume that the wind has been subject to a termination shock 
as to acquire a randomised direction of particle momenta (Gallant & Achterberg 1999). The 
number density of electron-positron pairs is uqj = PNs( r )^Ns/(47rec) (ncj is the Goldreich- 
Julian density; Goldreich & Julian (1969)) where e is the electron charge. 
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We adopt 7 W ~ 10 6 as the relativistic factor of most of the wind particles (see e.g. 
Gallant & Arons (1994)), and the ratio of ions to electron-positron pairs to be a = n-JriGj ~ 
10~ 3 ; ions carrying most of the energy of the wind. For future purposes, we write the ion 
density as 

niW = ««Gj(— f , (10) 
r 

where the index (3 describes the radial dependency. 

The QN ejecta will interact with the parent NS wind. The randomised wind particles 
will be boosted with one shock crossing by a factor of 2Fq N ~2x 10 4 before they leave the 
region 1 . The maximum number of particles that can be accelerated is limited by the QN 
fireball energy and can be written as (here the energy of the pairs will be negligible, since 
energetically they are at much lower gamma factors due to synchrotron losses) 

where m p is the proton mass, Z can vary from 1 to 26 (iron) and e is the fraction of QN 
energy transferred via the QN shock to the pulsar wind particles . If all ions are protons and 
approximately accelerated to the energies of the ankle (around 10 18 ' 8 eV), then up to 10 45 
particles would gain that energy. Heavier ions will correspondingly reach higher energies with 
the same total relativistic factor 2Fq N 7 w , e.g. iron will easily gain energies over 10 20 eV. 
Therefore in this model at higher energies there should be a natural increase in heavier 
ions. Since the above follows closely calculations already presented in Gallant & Achterberg 
(1999), the corresponding spectrum also in our case can be shown to be dN/dE ~ E~ 2 or 
even flatter given its dependency on Tq^. We would like to emphasize that further studies 
of our model are needed in order to make better predictions of the spectrum. 

Such an injection spectrum might agree with AGASA measurements but let us recall 
that in the high energy region we discuss here, currently there seems to be a disagreement 
between the AGASA ground array (Takeda et al. 1999) and the HiRes fluorescence detector 
(Abu-Zayyad et al. 2002) which seems consistent with the GZK-cutoff. Clearly there is a 
need for much larger experiments such as Auger, EUSO, and OWL, that can increase the 
number of detected events by one or two orders of magnitude before the injection spectrum 
is known conclusively. As we have said, for now, our model lacks the details to predict the 
exact spectrum. 



1 Particles do not have time to re-isotropise upstream before being overtaken by the QN shock; the 
acceleration consists of only one cycle. 
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The number of particles per unit volume from QNe in all the galaxies is J QN = 
(c/47r)iV wind n 9 T loss (i? a )z/Q N , where T\ oss {E a ) is the residence time of a particle at the an- 
kle initially injected at higher energy (Arons 2003). The QN rate per galaxy is given by z/qn 
while n g is the galaxy density. With z/ QN = 10~ 6 yr _1 and n g = 0.02 Mpc -3 , this implies 

jQN {E > 10 18.8 eV) ^ x 1Q _ 18 cm _ 2 gter _ 1 (l2) 

as compared to the observed value J ohs {E > 10 18 8 eV) ~ 3 x 10" 18 cm~ 2 s _1 ster -1 (1 event 
per square kilometer per year; Lawrence, Reid, & Watson (1991); Takeda et al. (1998)), 
assuming Z — 1. For particles at energies higher than 10 18 ' 8 eV we used an average Ti oss ~ 
1 Gyr (Biermann&Strittmatter 1987; Berezinsky & Grigoreva 1988; Bertone et al. 2002; see 
also Figure 5 in Sigl 2004b). 

The QN flux calculated above is short of the observed one unless the QN is very efficient 
in transferring energy via the QN shock into relativistic particles. However the presence of 
large scale magnetic field complicates the interpretation of UHECRs data and could lead 
to an overestimate of the observed flux. More specifically, the energy spectrum of particles 
emitted from a (nearby) source depends very strongly on the relative position of the observed 
source and the observer to the direction of the large scale field. In some cases the flux can 
be enhanced by a factor of ~ 100 for energies below 10 19 eV (see Stanev et al. 2003 and 
references therein). The observed flux and spectrum might not reflect the ones at the source 
making uncertain a direct comparison with the numbers derived above. 



4-1.1. Acceleration timescale 

The ejecta remain in expansion phase until the shock starts to decelerate. The transition 
between these two phases occurs at a radius -Rqn.cI at which the energy in the swept-up 
material becomes of the order of the energy released in the QN, i.e when a large fraction 
of the fireball energy has been used to reaccelerate the wind particles (Blandford & McKee 
(1976); Gallant & Achterberg (1999)). Using Eq.(10), 

i? QNd * ( -^qn(3 — (3) 

\ 47ran G j-RQ N 2r 2 7 w Zm p c 2 

which gives -Rqn.cI ~ 5.5 x 10 13 /\^Z cm for (3=1. If the ion density decreases as ~ 1/r 2 
{(3 = 2), the fireball will be in the expansion phase much longer and the distance should be 
defined by using interstellar matter density together with the wind (we estimate -Rqn )( i > 10 15 
cm for (3 = 2). These numbers will define the timescale of the highest energy cosmic ray 
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component acceleration. The case with (3 ~ 1 would more closely follow the magnetic 
field value and therefore the charge density of the pairs. We use it as an example. The 
corresponding time is tQN,d = -Rqn,ci/ c ~ 1700/v^ s. If all wind particles are iron, then the 
timescale would be shorter, 330 s. One should notice however that the ion density profile 
in the pulsar wind bubble is not known and therefore the timescale of the expansion phase 
could be much larger (with (5 ~ 2 up to weeks). 

4.2. Acceleration of the propelled wind 

The second cosmic ray component occurs when rapidly rotating millisecond pulsars with 
large magnetic fields, B > 10 14 G, undergo a QN explosion. In this case the born QS is in the 
propeller regime (see §3). The material flung away by the propeller expands as a magneto- 
hydrodynamic wind to reach the speed of light (with Lorentz factors up to 1000 at the light 
cylinder for B > 10 14 G; see Fendt & Ouyed (2004) for detailed calculations) 2 . The MHD 
nature of the propeller ensures that acceleration occurs at large distances from the star where 
synchrotron losses are likely to be minimum. As demonstrated in Fendt & Ouyed (2004) 
the magnetic field above the Alfven surface is predominantly toroidal. Such a geometry 
will allow the particles to escape freely along the poloidal direction in the acceleration zone 
without being deflected by the magnetic field lines. Particles with r prop > Tqn can reach 
eventually the QN shock, as to reach Lorentz factors as high as 2Fq N x r prop ~2x 10 7 . It 
is roughly 10 16 eV for protons and 10 18 eV for iron in one shock crossing. Since an efficient 
propeller needs strong magnetic field and a short period, this component is created in young 
pulsars soon after the supernova explosion or in binaries where the neutron star spins up 
(which also provides extra mass to trigger the QN). 

If the propeller works with the extreme accretion rate of rh — 10 28 g s _1 , then the rate 
of particle number is iV prop . = m/(Zm p ) ~ 10 52 /Z s _1 . That is, after one crossing, where a 
particle is boosted by a factor 2Fq N , the energy in time unit used to accelerate the propeller 
wind is iV2r prop m p c 2 ; only the particles with r prop > Tq N will be accelerated by the QN 
shock. This indicates that particle acceleration consumes the energy of the shock wave in ~ 1 
ms, and the shock dies out 3 . We expect the propeller to function efficiently in its early stages 
with most particles acquiring r prop > Tqn- In reality this timescale is longer because the 



2 Shvartsman (1970) first suggested that relativistic particles can be formed at the propeller stage by a 
rapidly rotating magnetic field (see also Kundt (1990)). 

3 Some of the ejecta may fall back allowing for disk and later planet formation around the newly born QS 
(Keranen & Ouyed 2003). 



-10- 



propeller wind particles cannot propagate straightforwardly to the shock (due to magnetic 
field and turbulence), their energies must have spread out and the most energetic particles 
reach the shock first; nevertheless, it is short compared to the timescale of the expansion 
phase. It also means that a tiny amount of the bulk of the propeller wind will be accelerated 
to energies of the order of 10 15 eV and at most up to 10 18 eV for those few particles that 
managed to get one more kick in the process. Here again, given the many similarities to 
the calculations and approach presented in Achterberg et al. (2001), we expect the boosted 
propeller wind to acquire a power law spectrum, dN/dE ~ E~ s with s ~ 3.2 — 3.3. 

Using the Milky Way dimensions with a radius of 15 kpc and a scale height of ~ 
1 kpc, we obtain the number density of galactic cosmic rays for one QN to be of the order of 
n MW ~ 4 x 10" 18 cm~ 3 /Z (10 15 eV/E)(E QN /10 53 erg). We assume a leaky box model (Garcia- 
Munoz et al. 1987; Simpson, & Garcia-Munoz 1988) for the galactic cosmic rays diffusing out 
of the galaxy where we adopt the leakage timescale expressed as T d = 3 x 10 7 yr (E Gc v I "ZY 1 ^ 
(see e.g. Webber 1998; Biermann et al., 2001 and Biermann & Sigl, 2001). The estimated 
time-averaged flux in our model is, 

J 2 QN (£> 10 15 eV) ~ fnMwTdi^QN 

47T 

3 x 10" ,10 eV 4/3 z/qn , e Eq^ -2-1 -1 /i A \ 

— 777z — ) ;r r — zr* — c m s sr , (14) 

Z 2/3 V E ) 4 - 6 yr- lA 10 52 erg ; ' V ! 

which is similar to the observed value of J ohs (E > 10 15 eV) ~ 3 x 10~ 10 cm' 2 s -1 ster' 1 (Bird 
et al. (1995); Candia (2003); Haungs (2004)) as long as the hyper-Eddington accretion rates 
can be accepted and if the accelerated particles are mainly protons. Notice that below the 
knee J^ N is too low to account for the observed cosmic ray flux; the shock is not energetic 
enough to accelerate more particles as to extend the spectrum to lower energies. Therefore, 
the steeper spectrum above the knee provided by the QN shock cannot dominate at energies 
below the knee. More important, the flux in Eq.(14) is averaged over time, t 3> I/Vqn- Above 
the knee the cosmic ray diffusion time is shorter than the QN occurrence, T d < r QN = 1/Vqn- 
For example, at the knee Td ~ 3 x 10 5 years for protons and it is Td ~ 3 x 10 4 years at the 
ankle, while the QN occurrence time is Tqn ~ 10 6 years. Therefore the propelled wind of 
the QN model can account for the galactic cosmic ray flux locally, i.e. in the vicinity and 
shortly after the QN event. Over the whole galaxy and in timescales longer than millions of 
years, QNe can contribute on average a few percent of the galactic cosmic rays around the 
knee. Closer to the ankle the contribution is negligible due to the very short diffusion time 
T d . More accurate estimates of the leakage times at high energies and QNe occurrence would 
be needed to make the above conclusions firm. In particular, the QN rate carries substantial 
uncertainties; among them the difficulty of determining the critical density for the transition 
to quark matter, the burden of not knowing which equation of state better describes the 
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parent NS, and the lack of a complete theory that can describe the QCD phase diagram that 
would describe the path followed by the NS during its transition to a quark star. We note 
however that recent studies show that QN rates can be as high 10~ 4 yr -1 galaxy -1 (Yasutake, 
Hashimoto, & Eriguchi 2004). which is above our fiducial value of 10 -6 yr _1 galaxy -1 . 

Finally, since the propelled wind consists of the crust material of the parent NS, we 
expect this cosmic ray component to be rich in heavier ions, up to iron. Their rigidity (E/Z) 
means that these heavy nuclei would leak out of the galaxy at a lower rate than protons. It 
implies an increase of the heavy nuclei to hydrogen ratio in the chemical composition over 
time for those cosmic rays originating in QNe. 



4.3. Quark Nova as a relativistic particle injector 

In our model the propeller injects relativistic particles. As discussed above, some of this 
propelled material can be re-accelerated by the QN shock. The remaining bulk will populate 
interstellar space. These particles with Lorentz factors T ~ 1 — 1000 (or energies 1-1000 GeV 
for protons) could later be reaccelerated by galactic supernova shocks. For the number of 
injected particles a simple estimate gives AT in j ~ E TOt /(10 GeV) ~ 10 52 erg/10~ 2 erg ~ 10 54 
that are accelerated to GeV energies by the propeller. Here E TOt is the maximum rotational 
energy of the QS. Using Eq. 14 this implies J ~ 10 -3 cm _2 s _1 sr _1 , which is a few percent of 
the observed value (see e.g. Abe et al. (2003)). In this case the particle density in the galaxy 
can reach umw = 4 x 10~ 13 cm~ 3 from a single QN using up the QS rotational energy. The 
spectrum of the population of these injected particles will be shaped by later encounters and 
acceleration by galactic supernova shocks. We note here as well that at most a few percent 
of the cosmic rays below the knee can be of the QN origin. 

The QN model predicts three cosmic ray components, two galactic ones and an extra- 
galactic one. This should in principle lead into transitions or steps in the cosmic ray flux. 
The galactic components will be subject to many SN shocks that will smoothen the discon- 
tinuity around the knee. Note however that these signatures will be dwarfed by the much 
higher contribution from other galactic cosmic ray sources. Having shown that QNe are 
candidate sources of UHECRs above the ankle, we next discuss how observational features 
can be explained within the QN model. 
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5. UHECRs from QNe in future cosmic ray detectors 

With a rate of 10~ 6 per galaxy per year and a galaxy density n g ~ 0.02 Mpc~ 3 one 
should expect about 0.08 QNe a year within a 100 Mpc sphere. This 100 Mpc region 
contains galactic and intergalactic magnetic field regions to which GZK-energy cosmic rays 
are subject to. This may lead to arrival time delay and therefore clustering of events around 
the source. The arrival timescales can spread from years to millions of years depending on 
the strength and the configuration of the magnetic field as well as the distance to the QN. 
Here we refer the interested reader to Sigl (2004a) for a recent discussion on the effect of 
the magnetic field on the ultra-high energy cosmic rays. Given these numbers, there may be 
a tiny energy window with the particles clustering within a few degrees towards the source 
(e.g, Kronberg 1994a and 1994b). This wishful possibility lead us to consider such events in 
the future detectors what we refer to as observational particle astronomy. Naturally, such 
observations will be better guided if the QN can be detected by other means as well (e.g. 
X-rays, gamma-rays). 

The energy lost by the UHECRs as they propagate and interact with the cosmic mi- 
crowave background is transformed by cascading into secondary GeV-TeV photons (Protheroe&Stanev 
1996). This TeV-gamma-UHECR trace could in principle be detected in the future (Cate- 
nese&Weekes 1999) and be used to test our model (see discussion in Akerlof et al. 2003 and 
references therein). However photons and UHECRs may have very different arrival times 
that are not easily quantified. This again calls for better understanding of magnetic field 
effects before such connections between the TeV photons and the related UHECRs can be 
made. Also, since protons can be accelerated up to 10 21 eV in our model, significant neutrino 
fluxes (with energies above 10 18 eV) can be generated (e.g., Waxman&Bahcall 1999; Engel 
et al. 2001). This is below the currently advertised threshold of 5 x 10 19 eV for EUSO and 
OWL and most of the potential events will go undetected (Halzen&Hooper 2002). Never- 
theless future neutrino detectors should be able to signal any neutrino-UHECR trace with 
the arrival direction of Ultra-High- Energy neutrinos as a good indication of the QN location. 
The above mentioned traces are the subject of another study of QNe as sources of UHECRs. 



5.1. Multiple events and clustering of UHECRs 

Assuming that the cosmic rays above the ankle are accelerated in QNe, we obtain 
roughly 10 44 particles above 10 19 eV per QN. The integrated flux F per unit area from one 
QN is 

F ^ 10 -io ( 100Mpc )2cm _ 2) (15) 
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where the Dqn is the QN distance. This is the total flux integrated over time including 
the spread in arrival time induced by magnetic fields, which can vary from years to millions 
of years (Sigl 2004a). In a detector like Auger this implies (integrated over the year and 
aperture and assuming that the QN is within the sky coverage) few hundred events per one 
QN at a distance of 100 Mpc. In a detector like EUSO and OWL one can expect one to two 
orders of magnitude more events per QN. If the time delay and thus the spread in arrival 
time is e.g. 1000 years, one could expect doublets and triplets around a given QN within 10 
year time in Auger, and equivalently dozens of clustered events with EUSO and OWL. For 
completeness, if the average time delay of cosmic rays above 10 19 eV is around 1000 years, 
in our model it would mean that UHECRs from roughly 100 QNe can presently be detected 
within 100 Mpc. 



5.2. Isotropy of UHECRs 

Unlike the arrival directions of UHECRs, the galaxy distribution is not isotropic within 
100 Mpc distance. This may mean that UHECR accelerators are not all within galaxies. 
It is natural to attribute UHECR sources to galaxies and the almost perfect isotropy, sofar 
observed both below and above the ankle (e.g. Torres&Anchordoqui 2004), to magnetic 
fields. Despite the uncertainties on the magnetic field strength it has indeed been shown 
that if our Local Supercluster contains a large scale magnetic field it can provide sufficient 
bending to the cosmic ray trajectories (Sigl, Lemoine & Biermann 1999; Farrar&Piran 2000). 
Heavy ions like iron, if present in the pulsar wind, will have unique consequences. Iron has a 
higher cutoff energy (Puget, Stecker, & Bredekamp (1976); Stecker, & Salomon (1999)) than 
lighter nuclei thus allowing a longer path of propagation in a tiny energy window allowing 
for more distance sources to contribute to the flux. Furthermore, being bent in galactic 
and intergalactic fields, heavy nuclei would acquire a more isotropic distribution of arrival 
directions than protons. Finally, heavy nuclei will naturally acquire higher energies than 
protons and contribute to the extreme end of the UHECR flux in our model. 

Of interest to the QN model, old population pulsars with large peculiar velocities 4 
could cover tens of Mpc distances from their origins in a Hubble time. Some of these run- 
away pulsars may undergo a QN explosion as a result of an increase of their core densities 
following spin down or accretion from the surrounding space. This would lean towards a more 



4 It is a well known fact that pulsars in our galaxy have velocities much in excess of those ordinary stars 
(Harrison, Lyne & Anderson 1993). It is reported that their transverse speeds range from to <~ 1500 km 
s _1 and their mean three-dimensional speed is 450 ± 90 km s _1 (Lyne & Lorimer, 1994). 
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isotropic distribution - the extent of isotropy in the runaway pulsar population remains to be 
determined - of the arrival directions of the UHECRs (the QN-shock boosted pulsar wind). 
This favours weak magnetic field (B < 10 10 G) millisecond pulsars in our model, with the 
assumption that these old isolated pulsars still sustain their wind bubbles. We note that 
these candidates are probably not born in the propeller regime following the QN explosion, 
given their weak magnetic field and long periods (see § 3). 

5.3. Injection spectrum and time-dependency 

For a given distant extragalactic QN source the more energetic particles should arrive 
first since the associated time delay induced by the intergalactic magnetic fields is short. 
Therefore it would appear as if there is a monoenergetic flux of particles from the source at 
one given time. Later, the lower energy particles will arrive with larger time delay (G. Sigl 
- private communication). 

In the case of galactic cosmic rays, the propelled wind consists of the crust material 
of the parent NS; the rigidity (E/Z) implies that the heavy nuclei would leak out of the 
galaxy at a lower rate than that of protons. We thus expect an increase of the heavy-nuclei- 
to-hydrogen ratio in the chemical composition over time scales of some ten thousand years 
for cosmic rays originating in QNe. Note, however, that this is effect is local and could 
have importance only if we detect a nearby quark star. If the star has undergone a QN 
phase in the past, the local cosmic ray composition may reflect the time passed since the 
QN explosion. 

6. Summary and conclusion 

The cosmic ray acceleration in the QN model consists of three different components as 
illustrated in figure 1. The first component is due to the acceleration of the high energy 
particles in the pulsar wind bubble by the QN shock. We predict few hundred UHECR 
events above 10 19 eV for the Pierre Auger detector per a distant QN, while some thousands 
are predicted for the proposed EUSO and OWL detectors. Magnetic fields can lead into 
clustering of the predicted events with timescales spread from years to millions of years. 

The second component stems from QNe with the compact remnant born in the pro- 
peller regime (pulsars with high magnetic field and small period undergoing a QN explosion 
immediately following a SN) ejecting relativistic particles with T ~ 1 — 1000. The particles 
with r prop > Tqn will eventually interact with the QN shock, as to reach roughly 10 16 eV 
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(protons) and 10 18 eV (iron) in one shock crossing. Given the energy of the QN shock, only a 
tiny amount of the propelled wind particles can be accelerated, i.e. a maximum of 10 50 par- 
ticles around the knee energy (see Eq. 14). The propelled particles that were not accelerated 
by the QN shock are injected into the galactic space and will eventually become accelerated 
by the supernova shocks in the galaxy. QNe as possible UHECR sources seem to account for 
the observed extragalactic flux and can contribute partially (a few percent) to the galactic 
cosmic rays. We conclude by stating that despite the fact that there is no guarantee that 
QNe occur in nature, the model possesses features that can be tested in future cosmic ray 
detectors. 
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Fig. 1. — Cosmic ray components in our model. The solid line illustrates the observed cosmic 
ray flux while the dotted line illustrates the contributions from the QN model. Above the 
ankle are the pulsar wind particles accelerated by the QN shock. Below the ankle and above 
the knee are the particles first propelled by the QN compact remnant and then boosted by 
the QN shock. Only the most relativistic propelled particles will interact with the shock 
which gets shut-off after a few milliseconds. Following this latter phase the QN compact 
remnant acts as an injector of relativistic particles which later would get accelerated by 
galactic SN shocks to form the component with energies below the knee. Note that QNe 
as possible UHECR sources seem to account for the observed extragalactic flux (above the 
ankle) and can contribute partially (a few percent) to the galactic cosmic rays (below and 
around the knee). We note however that the observed extragalactic flux (above the ankle) 
awaits experiments such as Auger, EUSO, and OWL to resolve current disagreement between 
AGASA and HiRes. 



